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Abstract We introduce a quantum dot orbital tight-binding non-equilibrium 
Green's function approach for the simulation of novel solar cell devices where 
both absorption and conduction are mediated by quantum dot states. By the 
use of basis states localized on the quantum dots, the computational real 
space mesh of the Green's function is coarse-grained from atomic resolution 
to the quantum dot spacing, which enables the simulation of extended devices 
consisting of many quantum dot layers. 
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1 Introduction 

Extended quantum dot superlattices (QDSL) as found in third generation so- 
lar cell devices [TJ|2] consist of a large number of weakly coupled quantum 
dots, with a number of involved atoms exceeding the limit of what may still 
be handled by even today's largest supercomputers. For that reason, atom- 
istic [3] and even microscopic continuum [4 simulations of QDSL absorbers 
make use of the symmetry and periodicity properties present in the ideal- 
ized structure. However, in real devices, the finite number of periods, large 
built-in fields and any kind of disorder make the electronic structure deviate 
considerably from the minibands found for the latter and hence disqualify such 
simplifications. Furthermore, even for narrow electronic minibands and gaps 
in the phonon spectrum, the charge carriers will be subject to various kinds 
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of scattering, causing dissipation and preventing purely coherent transport in 
the device. For the determination of the device characteristics it is thus de- 
sirable to derive a computational model that on the one hand can still be 
handled numerically by making use of the localized nature of the QD wave 
functions, and which on the other hand is able to cope with all of the above 
mentioned deviations from the idealized structure, thanks to being based on 
a quantum-kinetic theory of photogeneration, recombination and transport, 
including injection and extraction of carriers at contacts, similar to the non- 
equilibrium Green's function (NEGF) approach developed for quantum well 
solar cells in [SUM?]. The present paper introduces such a model and applies 
it to the investigation of (photo)carrier localization and extraction in selec- 
tively contacted silicon quantum dot absorber arrays under varying internal 
field conditions. 



2 Theoretical approach 

The NEGF theory that forms the core of the model provides a quantum version 
of the Boltzmann transport equation for the charge carriers [Sl CTTPiriT] , which 
includes tunneling transport and non-local rate terms for scattering, genera- 
tion and recombination, which are all essential processes in the photovoltaic 
regime of device operation. For the representation of the electronic system, the 
present hybridization approach uses the perturbative expansion of the QDSL 
wave function in terms of the eigenstates of the isolated dots. The resulting 
molecular orbital approach can be interpreted as a tight-binding (TB) theory 
with quantum dot orbitals replacing the atomic orbitals, and is in spirit similar 
to the first NEGF models of quantum well superlattices for quantum cascade 
laser simulations [P21IT3] and to NEGF models of polaron transport in QD 
superlattices [2]. The orbitals m = 1, 2, . . . , M of QD i = 1,2,..., Nd with 
energies £i m are used to represent the field operators for single particle quan- 
tum dot states, ^(r,i) = J^. ipi n (r)di n (t), where dm annihilates an electron 
in state \n) on the QD at position Rj. The non-interacting nearest-neighbor 
tight-binding Hamiltonian in this QD orbital basis is 

M N D M 

(*). (1) 

m,n— 1 i—1 m—1 

where are nearest-neighbor sites, t is the hopping matrix, n the den- 

sity operator and £i m = Si m + Ui, with fZj the average value for the Hartree 
potential of Coulomb interaction at the dot position. The heterostructure po- 
tential does no longer appear explicitly, since it has been considered in the 
determination of the TB-parameters e and t, 

Um tj n = ~ [ d 3 r^ m (r)AU(r)^ n (r) (2) 
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with AU(r) = Hsl{y) — J^Ri Hqd{y — R-i), where Hsl is the superlattice 
Hamiltonian and Hqd is the Hamiltonian of the isolated QD. The Hamilto- 
nian ([I]) is used in the equations for the steady state non-equilibrium Green's 
functions for the device states, which read 



-l 



G R (E) = {G R (E)}~ - S RI (E) - S RC {E) , (3) 

G$(E) = G R (E)[S^ I (E) + E^ C (E)]G A (E), (4) 
G R (E) = \{E + ir,)l - Ho]" 1 , G A (E) = [G R (E)}\ (5) 

where G = [G lmJn ] with G im j n (E) = J dTe tET/h G imJn (T), r = t' - t, and 
Gi m j n (t,t') — — ^(Tcldimitydjnit')}) the quantum-statistical non-equilibrium 
average defined on the Keldysh contour C [5]. In the above equations, the 
self-energy terms S aI (a e {sg, R, A}) describe the renormalization of the 
single charge carrier Green's function due to the interactions relevant for the 
photovoltaic device operation, i.e., coupling to photons, phonons and other 
charge carriers, and are obtained within many-body perturbation theory on 
the level of the self-consistent first Born approximation. As these self-energies 
are functionals of the carrier Green's functions, they need to be computed 
self-consistently with the latter. The additional self-energy S aC encodes the 
coupling of the QD states to the states of the contacts, which in this case are 
assumed as bulk electrodes. Details on the general derivation and explicit form 
of the self-energies can be found in [7] . 

Once the NEGF have been determined, they immediately provide any phys- 
ical observable on the single particle level. For instance, the density of states 
is obtained from the NEGF and the QD eigenf unctions via 

9{t,E) =53E^«WS»WV'i«(r). (6) 

i,j m,n 

where A = i(G R — G A ) is the charge carrier spectral function. In terms of 
the tight-binding intcrdot hopping and Green's function elements, the current 
flowing between two quantum dots can be expressed atQ 

Ji = ~ J dE tr {t« +1 G< U (E) - G< +1 (E)t l+ll } . (7) 



3 Implementation for Si-SiC>2-SiC QD heterostructure absorbers 

Since the wave functions and energies of the quantum dot eigenstates are de- 
termined separately for each individual dot, the corresponding computational 
domain is small enough to allow, in principle, for the use of accurate ab-initio 
methods. In this paper, however, since the main focus is on the general formal- 
ism, the orbitals are approximated via the superposition of three separate ID 



1 The trace is over orbital indices. 
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Fig. 1 (Color online) (a) Schematic structure and (b) resulting band diagram of ap-i-n- QD- 
array absorber solar cell device. The resulting net potential drop Vw in the absorbing region 
is determined by the combination of the built-in potential due to doping and the terminal 
bias voltage Vf,i as , the latter corresponding to the separation of the chemical potentials fin,p 
at the contacts. 



solutions for a simple one-band effective mass model. In order to agree with 
full real space calculations, the local variation of the electrostatic mean-field 
potential in principle needs to be included in the solution of the Schrodingcr 
equation, and to be computed self-consistently with the carrier density via 
Poisson's equation. In this first approach, a linear potential drop resulting in 
a constant field is assumed and the orbitals are computed only once for the 
flat-band case. This means that quantum-confined Stark effects resulting from 
the distortion of the orbitals due to the local field are not considered at the 
present stage. 

The structure under investigation is shown schematically in Fig. [lja). It 
consists of decoupled arrays of 20 cubic silicon quantum dots of 3 ran side 
length embedded laterally in SiC>2 and vertically separated by thin SiC layers 
of 1 nm thicknes^] The QD array is coupled to contacts on both ends, and the 
regions adjacent to the contacts are doped to provide a built-in field. The non- 
selfconsistent band diagram is shown in Fig.[l|b). The internal field depends on 
the built-in potential Vm, determined by the doping, and the bias voltage Vbias 
applied at the right contact, the latter corresponding to a finite separation of 
the chemical potentials associated with the equilibrated electrodes. In order to 
prevent leakage current flow and to enable charge separation even in absence 
of an internal electric field, minority carrier contacts are closed artificially, the 
contacts thus becoming carrier-selective. 

At this stage of implementation, the optical modes considered are those 
of a homogeneous medium and with monochromatic occupation, and only 



2 A similar structure was recently fabricated as a candidate for an all-silicon tandem solar 
cell component |15j . 
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radiative recombination is considered. Dissipation of energy is modelled via 
coupling to a single optical (bulk) phonon mode described within the defor- 
mation potential approximation. The energy of the optical branch was chosen 
at E p hon = 0.06 eV. Additional broadening of the electronic states is obtained 
from elastic scattering with acoustic phonons, using again the appropriate bulk 
deformation potentials. 



4 Numerical results 

Fig. [2] shows the local as well as the spatially integrated density of states 
(LDOS and DOS) of the selectively contacted QD-array structure from Fig. [I] 
for vanishing internal field {Vu — Vbias = 0), in the case where only the lowest 
QD orbital is considered. Due to the finite number of dots, the density of 
states shows discrete maxima, but as a consequence of the flat band situation, 
the states are completely delocalized over the structure, such that on every 
QD site, the multiplicity of states amounts to the number of QDs, as can be 
recognized in the magnification of the shaded region (Fig. [3]). The difference 
between electrons and holes is due to the differing values for effective masses 
and barrier height. While the weak mutual coupling of QD states results in 
a very narrow "miniband" in the active device region, the LDOS is strongly 
broadened close to the carrier selective contacts, which is due to the efficient 
coupling to the broad-band bulk electrode DOS. This hybridization with bulk 
electrode states is also reflected in the integrated DOS. 

In the flat band situation, extraction of photogenerated charge carriers is 
efficient, since all the absorbing states are directly connected to the contact 
states. However, the situation changes drastically in the presence of a finite 
field, as displayed in Fig. |4j The width of the " miniband" increases in the in- 
tegrated DOS, but in the LDOS, the number of states per QD with significant 
spectral weight is reduced, since even for a field as weak as 25 kV/cm, which is 
close to the values encountered in real devices, the states are spatially localized 
on a few QDs. This carrier localization, clearly seen in the magnification shown 
in Fig. [5] prevents a direct ballistic extraction of photogenerated charge carri- 
ers and thus increases the significance of scattering assisted transport ranging 
from sequential tunneling to hopping between single dots. 

The restriction to the lowest orbital leads to an underestimation of the 
overall inter-dot coupling. In order to study the extraction of excess charge 
generated under illumination, the single-band hopping matrix element is in- 
creased to 0.05 eV. To investigate the transport of photogenerated charge 
carriers under conditions were the states are no longer completely delocalized, 
i.e., in the situation where inelastic scattering is crucial to enable carrier flow 
between quantum dots, the local photocurrent spectrum is computed for grad- 
ually increasing internal field. The result is displayed in Fig. [6] for monochro- 
matic illumination at a photon energy of 1.67 eV, corresponding roughly to 
the effective energy gap. Even though the current spectrum is strongly mod- 
ified by the internal field, the energy integrated current differs only slightly 
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Fig. 2 (Color online) Local density of states (LDOS) and integrated DOS of the selectively 
contacted QD-array structure from Fig. ^ for vanishing internal field {Vm = 0, flat "mini- 
band"). Only the contribution from the lowest orbital is shown. While the weak coupling 
of the QDs results in a very narrow band, the LDOS close to the carrier selective contacts 
shows a characteristic broadening due to hybridization with the states of the bulk electrode. 
The magnification of the shaded area is shown in Fig. [3] 




position z [nm] density of states [a.u.] 



Fig. 3 LDOS and integrated DOS at vanishing built-in field (shaded area from Fig. [2 i . 
In this situation, the states are maximally delocalized and extend over entire multi-QD 
structure. The discrete set of maxima in the LDOS on the QD positions reflect both coupling 
and finite number of QDs in the structure. 
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Fig. 4 (Color online) Same as Fig.|2j but for a finite internal field of 25 kV/cm. The field 
splits the quasi-miniband into a broader level structure with reduced local energetic overlap. 
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Fig. 5 (Color online) Detailed view of the situation at finite internal field of 25 kV/cm 
(shaded area in Fig. [5}. The field has a strong localizing effect on the QD-array states, 
which are spatially confined to a small number of neighboring QD. Without the presence 
of higher states, ballistic transport to the contact is no longer possible in this situation, 
i.e., energy dissipation via inelastic scattering processes is required to enable photocarrier 
extraction. 



for different field values, and the overall current, i.e. the sum of electron and 
hole contributions, is perfectly conserved over the entire device. The strong 
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Fig. 6 (Color online) Photocurrent spectrum at photon energy of 1.67 eV for zero applied 
bias and different values of the built-in field. Due to strong electron-phonon interaction, the 
photogenerated electron-hole pairs relax to the effective band edge within the derealization 
range of a few quantum dot spacings, and photocarrier extraction remains possible even 
in absence of ballistic current contributions. In spite of the large variation in the current 
spectrum for different internal fields, the energy integrated current differs only slightly, and 
the total current (electron+hole) is perfectly conserved over the entire device. 



electron-phonon scattering leads to a fast relaxation within the few QDs over 
which the states are delocalized, such that photocurrent flow follows the ef- 
fective band edge. Thus, efficient photocarrier extraction is recovered at the 
expense of the energy dissipated in the inelastic scattering process. 

5 Conclusions 

The QDTB-NEGF approach extends the applicability of powerful quantum- 
kinetic methods to optoelectronic devices based on electronically coupled quan- 
tum dot structures, with arbitrary internal fields and spatial arrangement. 
Application to an array of weakly coupled quantum dot reveals the significant 
carrier localization for finite electric fields in the absorbing region. For stronger 
coupling, photogenerated carriers are extracted efficiently also at considerable 
internal fields via phonon-assisted sequential tunneling and inter-dot hopping. 
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